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Abstract 

We investigate experimentally the drying of a small volume of perfectly wetting liquid on two crossed fibers. We 
characterize the drying dynamics for the three liquid morphologies that are encountered in this geometry: drop, 
column and a mixed morphology, in which a drop and a column coexist. For each morphology, we rationalize our 
findings with theoretical models that capture the drying kinetics. We find that the evaporation rate significantly 
depends upon the liquid morphology and that the drying of the liquid column is faster than the evaporation of 
the drop and the mixed morphology for a given liquid volume. Finally, we illustrate that shearing a network of 
fibers reduces the angle between them, changes the morphology toward the column state, and therefore enhances 
the drying rate of a volatile liquid deposited on it. 


1 Introduction 

Evaporation is ubiquitous in natural phenomena such as 
the transpiration in plants [T] , the perspiration of humans 
and animals, and the dynamics of sea spray [2]. Evap¬ 
oration also impacts many technological applications in¬ 
cluding aerosol drying, food dehydration, and coatings. 
However, modeling of the drying dynamics requires a de¬ 
scription of the geometry of the liquid, which is generally 
related to the wetting properties of the substrate. 

When a drop is in contact with a substrate, the evapo¬ 
ration dynamics strongly depends upon the wetting prop¬ 
erties. The simplest situation of a drying drop consists of 
isolated aerosol droplets or non-wetting liquids mi QIC] 
for which the evaporation rate is set by the radius of 
the droplet laE]. For wetting liquids, the effect of the 
contact line must be taken into account [iinilin]- In¬ 
deed, the contact line generates a stronger evaporation 
flux because of a “tip effect” [in d ini d Ell [n]. 
Furthermore, other properties of the substrate affect the 
drying such as the deformation of soft substrates m or 
drop confinement between two plates [18] . 

So far, most experiments address the drying of droplets 
resting on a plane but the drying of a drop on other sub¬ 
strates with different geometries remains largely unex¬ 
plored. For example, the evaporation of a liquid volume 
between fibers has industrial applications but remains 
poorly explored. Indeed, networks of fibers such as fab¬ 
ric, filters and paper are a collection of crossing fibers. 
In particular, at the node of two fibers, the liquid can 
adopt three different equilibrium morphologies that have 


been studied only recently [19]. 

Previous studies have investigated the morphology 
adopted by liquid on fibers. To understand the gen¬ 
eral situation of a complex three-dimensional network 
of fibers |2Q], the wetting has been studied in different 
fiber configurations. When a liquid is coated on a single 
fiber mi on 03], depending upon its rheological proper¬ 
ties mu QS] , the liquid can undergo a Rayleigh-Plateau 
instability, which results in a series of droplets [23]. A 
liquid droplet placed on a single fiber can adopt a bar¬ 
rel shape, a clam-shell shape or detach as gravitational 
forces dominate capillary forces miiiEiiis]- 

Other studies have investigated geometries involving a 
pair of fibers. Observations of the equilibrium morphol¬ 
ogy of a wetting liquid on two parallel fibers has revealed 
that two morphologies can be adopted by the liquid. 
When the distance between two fibers is large, the liquid 
is in a drop morphology in contact with both fibers. Be¬ 
low a critical distance between the two fibers, the fluid 
spreads into a column that spans the gap between the 
fibers mni on. Consequently, the interfacial area per 
unit volume of liquid increases significantly when the liq¬ 
uid changes morphology from a drop to a column, which 
speeds up the evaporation of the liquid [32] • More re¬ 
cently, the investigation of the more complex situation 
of crossed fibers has revealed a rich morphology diagram 
as illustrated in figure [TO]. Three morphologies can be 
obtained depending upon the volume of liquid and the 
tilt angle S between the fibers. 

In this paper, we investigate the drying of a completely 
wetting liquid on crossed touching fibers with different 
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oration speed can be defined as [32] 
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Figure 1: Pictures of the three possible morphologies for 
increasing tilt angle 5 of the fibers: column (red), mixed 
(blue) and drop (green). The fiber diameter is 2a and 
the total column length is 2 L. 


1 dm(t) 
Pi Sp dt 


( 1 ) 


where Sp = ttR^ is the surface area of the petri dish 
and the liquid density. The measured evaporation 
speeds are Ve = 3.9 x 10“® m/s and Ve = 3.6 x 10“® 
m/s for OMTS and DMTS, respectively. Note that it is 
consistent that the higher molecular weight (i.e. higher 
viscosity) liquid evaporated more slowly. 


tilt angles. We report experimental results for all mor¬ 
phologies adopted by the liquid on crossed fibers, i.e. the 
drop, the column and the mixed morphology. For each 
of the morphologies, we provide a model to describe the 
time evolution of the liquid shape during drying. We 
demonstrate that the tilt angle controls the efficiency of 
liquid evaporation and apply our findings to the evapora¬ 
tion of liquid on networks of fibers. The observation that 
the tilt angle between fibers influences the evaporation 
rate indicates that mechanical effects can affect drying 
of wet fibrous materials. 


2 Experimental section 

Two nylon fibers of circular cross-section and of the same 
diameter 2a (0.28 or 0.30 mm) are clamped on two hor¬ 
izontal rotating microcontrollers (PROl, Thorlabs). The 
tilt angle 5 (see figure between the two fibers is ad¬ 
justed using the microcontrollers. We deposit a liquid 
drop of volume V € [0.5; 4] pL using a micropipette on 
the node of the two fibers. Pictures are taken using Nikon 
cameras (D5100 and D7100) and 105 mm macro objec¬ 
tives to observe both the side and the top views. The 
tilt angle S between the two fibers is measured directly 
from the top views by image processing with an uncer¬ 
tainty of 0.02°. To ensure that the drying conditions are 
controlled, the entire setup is placed in an acrylic box to 
limit air convection. 

For all of the experiments presented in this pa¬ 
per, we use completely wetting volatile silicon oils 
purchased from Sigma-Aldrich: octamethyltrisiloxane 
(OMTS) with a kinematic viscosity 1 mm^/s, a density 
Pi = 818 kg/m^ and a surface tension 7 = 17.4 mN/m, 
decamethyltetrasiloxane (DMTS) with a kinematic vis¬ 
cosity 1.5 mm^/s, a density pt = 854 kg/m^ and a surface 
tension 7 = 18.0 mN/m. These silicon oils are composed 
of a monodisperse oligomer that ensures a single drying 
kinetics. 

The evaporation speed Ve is obtained by recording 
the mass evolution m(t) of a layer of liquid placed in a 
petri dish of radius Rp = 1.6 cm and weighed using a 
high-precision scale (XS105, Mettler Toledo). The evap- 


3 Results 

3.1 Observations 

The drying of liquid lying on crossed fibers, involves three 
different morphologies (figure[^ [T3] . Indeed, as shown in 
figure the morphology adopted by the liquid depends 
upon the tilt angle <5 and the drop volume V (see supple¬ 
mentary information video) and changes in time during 
drying. In particular, as evaporation proceeds, the vol¬ 
ume decreases and the new equilibrium morphology can 
be read on a morphology diagram for a given value of 
the tilt angle (see examples in figure [^b)). Therefore, 
we observe that liquids initially in the column morphol¬ 
ogy remain a column throughout the entire evaporation 
process while liquids initially in the mixed morphology 
switch to the column state. When the initial morphol¬ 
ogy of the liquid is a drop, the liquid can be assumed to 
remain a drop since the transition to the column state 
occurs at a very low volume. These states and transitions 
are summarized in figure [^b). 

3.2 Theoretical background 

Denoting py the vapor concentration in the gas surround¬ 
ing the liquid volume, Tick’s second law is 

^ = Dy^V^Py, ( 2 ) 

where £>m is the diffusion coefficient of the vapor in the 
gas phase. Considering a characteristic lengthscale £ and 
an evaporation speed Ve, the Peclet number Pe |33| is 
the ratio of the diffusive time of the vapor in the gas 
phase te = ^jDm to the evaporation time te = IjVE'- 


The diffusion coefficient Dm is derived from the kinetic 
theory of gases. Dm ~ 10“® m/s^ [34] . 

In our experiments, the characteristic length scale i is 
the typical size of a drop, which is also the order of the 
capillary length = \/^Hpig) — 1-5 mm, where g is 
the gravitational acceleration. The characteristic Peclet 
number is thus small, Pe <C 1: the gradient of vapor in 
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the gas surrounding the liquid is established much faster 
than the evaporation of the liquid. Consequently, the 
time derivative in equation (|^ is negligible. Therefore, 
the drying dynamics can be described by the Laplace 
equation: 


= 0. (4) 

We define the flux jo as the volume of liquid that evap¬ 
orates per unit area and unit time, and write [5] 

d Py 
dr 

where tq denotes the position of the interface. Assuming 
that pv(r = ro) = Psat, the vapor density at saturation, 
equation © provides the typical magnitude for the dif¬ 
fusive flux jo 


PUo 

Dm 


(5) 


jo - ^ (6) 

where D = Dm Psat jPi and £ is the length scale of the 
vapor concentration gradient. 

In the present paper, we assume that the flux jo is 
uniform over the air/liquid interface i.e. we neglect the 
increase of the evaporation flux near contact lines. Thus, 
we can write that 


dl^ 

dt 


- -jo^ 


DA 


(7) 


where V is the volume of liquid and A the total interfacial 
area. Equation Q shows that the drying rate depends 
upon the morphology of the liquid through the interfacial 
area A and the lengthscale of the vapor concentration 
gradient £, which is set by the liquid geometry. Next, we 
analyze different liquid morphologies and evaluate these 
two parameters, A and £, to model the drying dynamics. 


3.3 Drying of a droplet at the node of 
two crossed fibers 


Figure 2: (a) Time evolution of a drying liquid at the 
node of two fibers for three different tilt angles S. The 
initial volume is 2 pi (OMTS) and the fiber diameter is 
2a = 0.30 mm. Black arrows indicate the position of 
the liquid column extremities. Colored lines on the left 
of each picture refer to the liquid morphology (see figure 
[^. (b) Phase diagram of the liquid morphology as a 
function of the tilt angle S and the dimensionless volume 
V = VjcA. Arrows indicate the evolution of the system 
in the phase diagram during the drying process for the 
three tilt angles. 


The drying dynamics of a spherical drop are set by the 
radius of the drop i?, i.e. £ = i? [5]. Thus, the flux is 
jo — DjR and the drop radius decreases according to 
R{t) = y^2D{tf — t) with tf being the total evaporation 
time. 

For a drop on a single fiber in an axisymmetric barrel 
shape [27l |28] , the drop is characterized by its length L 
and its height h (figure [^a)). The height h is defined 
as the distance between the drop apex and the fiber sur¬ 
face. Assuming that the typical lengthscale of the vapor 
concentration gradient is £/2 [32], the evolution of the 
drop height is described by [35]: 


dt L{dV/dh)' ^ ' 

The interfacial area A and the volume V can be obtained 
analytically from Carroll’s derivations [S^ [35]. 
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Figure 3: Definition of the fluid height, h, (a) for a sin¬ 
gle fiber and (b) for crossed fibers, h = {hi + /i 2 )/ 2 . (b) 
Drying dynamics of drops in three configurations: single 
fiber and crossed fibers with tilt angles <5 = 65° and 90°. 
The plot represents the evolution of the height h/hmax 
with {tf — where hmax is the initial height and t/ is 
the total duration of the drying process determined ex¬ 
perimentally. In each case, two liquids are used (OMTS 
and DMTS). The initial volume is 1 and the fiber di¬ 
ameter is 2 a = 0.30 mm. The black line is the theoretical 
prediction given by equation ([^. 





Figure 4: Drying of a long column on two fibers with 
a tilt angle S = 0.45°. The initial length of the liquid 
(DMTS) column is Lq = 7 cm and the fiber diameter 
is 2a = 0.30 mm. The inset shows the raw data, the 
length of the liquid column L as a function of time tf — t. 
The main plot represents the data with ln(L/Lo) versus 
\/{tf —t)/tf. The black line is the theoretical predic¬ 
tion (equation ([^) for the drying of long liquid columns 
between parallel fibers |32) . 

In the case of crossed fibers, the definition of the fluid 
height needs to be adapted as illustrated in figure [^b) . 
The distance is equal to the average along the two bi¬ 
sectors. It is noteworthy that the lengthscales of the 
droplet on single and crossed fibers are similar. In figure 
[^c), we report the evolution of the drop height h rescaled 
by its initial value hmax as a function of {tf — t^jtf. 
The total duration of the drying process tf is determined 
experimentally from the direct observation that no liquid 
remains on the fibers. The numerical solution of equation 
([^ is calculated with an iterative method as explained in 
the appendix of reference [32] . In both single and crossed 
fibers configurations, the drying dynamics follow a trend 
that is in a good agreement with the numerical solution. 

In addition, we observe good agreement between ex¬ 
periment and theory for the two liquids with different 
volatilities in each configuration as reported in figure 
Thus, the results validate that the vapor concentration 
gradient, and thus the dynamics, is set by the size of the 
droplet. Consequently, the evaporation of a drop at the 
node of two crossed fibers can be fairly well described 
by a model similar to the one developed for a drop on a 
single fiber. 

3.4 Drying of liquid columns 

We now describe the drying of the liquid in the column 
morphology, which is observed for small enough tilt an¬ 
gles or volumes of liquid (figure]^. In such situations, 
the Bond number Bo = piga^/j that compares gravi¬ 
tational forces and surface tension effects is of the order 
10“^. Thus surface tension effects dominate and we shall 


4 














Figure 5: Schematic of the cross section of a liquid in a 
column morphology between two crossed fibers. 


assume in the following that the liquid maintains a sym¬ 
metrical shape despite gravitational forces. 

3.4.1 Small tilt angles and long columns 

We first consider the case of very long columns obtained 
for small angles 6 (figure[^ . We define the length L of the 
column from the node of the fibers to the liquid column 
extremity as illustrated in the figure [l] We represent the 
time evolution of the length L obtained for very small tilt 
angle, 6 = 0.45°, and a column of initial length Lq = 7 
cm, in figure]^ Note that in the hnal stages of the drying 
process, the length of the column cannot be measured 
near the node. 

In order to model the drying, we assume that the dry¬ 
ing of liquid columns between parallel [32] and crossed 
fibers are similar. We assume that the liquid interface 
has an ellipsoidal shape with a large aspect ratio, i.e. 
b L where b = Ltan{S) is the distance between the 
fibers. Thus, we set the lengthscale of the vapor concen¬ 
tration gradient C « bln{L/b). As derived in [52], we 
can show that 


ln(L/Lo) oc (9) 

Figure [^ shows the agreement between the experimen¬ 
tal observations and the theoretical prediction given by 
equation However, in practice, small tilt angles are 
relatively rare. Indeed, in a randomly oriented network 
of hbers, the angle distribution can be assumed to be 
uniform. It is therefore important to consider columns 
formed at relatively larger angles. Indeed, as the angle 
5 increases, the column length L decreases rapidly [19] 
and the drying dynamics changes. 

3.4.2 Modest tilt angles and short columns 

For larger tilt angles i5, the aspect ratio of the liquid 
column can no longer be considered large. In this case, 
the time-evolution of the volume of liquid in the column 
state can be expressed as 

nL{t) 

V{t) = 2 / A{y)dy, (10) 

Jo 


where the factor 2 corresponds to the two half-columns 
on each side of the node, and A{y) is the cross-sectional 
area of the column (figure [^. Denoting s{y) the arc 
length of the liquid cross-section, the total interfacial 
area is given by 


fL{t) i-Lit) 

A(t) = d / s(y)dy = 4 / {-k - 2 a{y)) R{y) dy, 

JO Jo 

( 11 ) 

where a{y) and the radius of curvature Riy) of the liquid 
are defined in figure]^ and detailed formulae are given in 
[12]. Substituting equations (10) and (11) into equation 
([^, we obtain 

\ 2]j rL{t) 

^ A(y)dyj=--^y {ir - 2 a{y)) R{y) dy, 

( 12 ) 

which simplifies to 



2D 


dL{t) _ 

dt L{t) A [L{t)] 


rm 


(tt -2a{y)) R{y)dy = 0. 

(13) 

Experimentally, we choose conditions for which a col¬ 
umn is observed at angles 6 G [5,16]° and we measure the 
time evolution of the length of the column for different 
initial volumes. These results are rescaled with the initial 
column length and total evaporation time and reported 
in figure]^ We note that the data all collapse on a master 
curve and are in good agreement with the numerical so¬ 


lution of equation (13). In summary, we have succeeded 


in describing the drying behavior of long columns (for 
small tilt angles) and short columns (for larger tilt an¬ 
gles) successfully with our model although the transition 
between these two regimes may be more complex to cap¬ 
ture. However, for practical application in an array of 
fibers, the tilt angle will remain large enough to apply 
the second situation or small liquid column. 


3.5 Drying of liquid in the mixed mor¬ 
phology 

For intermediate tilt angles between fibers, we observe 
a mixed morphology (figure]^, which consists of a drop 
on one side of the node and a liquid column on the other 
side [19] . In figure [^a) , we show the time evolution of 
the liquid length on both sides. The dynamics can be 
split into two regimes. Initially, the liquid adopts the 
mixed morphology state with a large drop in equilibrium 
with a small column. During the drying process, the 
volume of the drop decreases while the column elongates 
slightly to balance the Laplace pressure [12]. This first 
regime ends when the two lengths are comparable. In the 
second regime, the drying of two liquid columns (figure 
[^b)) follows the dynamics described in section 

In the first regime, the drop height is measured from 
the side view. The shape is slightly asymmetric due to 
gravity effects since the Bond number Bo = pegV^^^j^ ^ 


3.4.2 
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Figure 6: Evolution of rescaled length L/Lq for the col¬ 
umn state on crossed fibers with rescaled time (tf — 
for various initial volumes of OMTS. The fiber diameter 
is 2 a = 0.28 mm and the tilt angle is 5 = 10°. The solid 
line is the time evolution of the length as predicted by 
the numerical solution of equation (13). 


10“^ is not very small, in these experiments. Neverthe¬ 
less, as shown in figure [^c), the height evolution be¬ 
haves similarly as drops on single fibers or at the node of 
crossed hbers (figure]^. Consequently, we show that the 
evaporation rate in the mixed morphology state is dom¬ 
inated by the drying of the drop. When the liquid vol¬ 
ume reaches a critical value, the liquid adopts a column 
shape. Figure [^b) illustrates that the time evolution of 
the column length is well described by the equation (13) 
for short columns, which is consistent with that of small 
columns described in the previous section. 
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Figure 7: Drying dynamics starting from a mixed mor¬ 
phology. (a) Evolution of the liquid length on the column 
and the drop sides with rescaled time (tf — where 
tf is the total drying time. The length is normalized by 
Lq, the length of the columns at the morphology transi¬ 
tion. The red points represent the drop side and the blue 
points the column side, (b) Time evolution of the col¬ 
umn lengths in a situation similar to figure (c) Time 
evolution of the drop height as in figure [^b). The fiber 
diameter is 2a = 0.30 mm and the tilt angle is S = 20°. 
Times and are the drying times for the column and 
the mixed morphology states, respectively. An initial 
volume of 2 fii of DMTS is deposited at the node. 


3.6 Network of fibers 

Building on our understanding of the evaporation dy¬ 
namics in the three morphologies, we compare drying 
measurements for different tilt angles. In figure we 
present the time evolution of the liquid length for the 
three different initial morphologies. The data illustrate 
the transition between the liquid morphologies as well as 
their respective durations of drying. We also note that 
the drying kinetic for liquid columns (blue points in fig¬ 
ure]^ follow the same behavior, independent of the ini¬ 
tial morphology. As discussed previously, due to both the 
surface area and the characteristic vapor concentration 
gradient £, the evaporation rate is larger for a column 
than it is for a drop. 

Now, we apply this finding to the drying of liquid at 
the nodes of a network of fibers. To directly observe the 
effect of the fiber tilt angles on the drying rate, we make 
a 2D network of crossed fibers. The network is com¬ 
posed of 16 nodes and each node is about 1 cm from the 
neighboring nodes, so that the nodes may be considered 



Figure 8: Starting from the three possible morpholo¬ 
gies, the evolution of the length L/Lq is represented as a 
function of {tf — The shapes of the symbols corre¬ 

spond to the initial morphology and the color represents 
the morphology for the given point. Data for liquid ini¬ 
tially in columns (-f) is replotted from figure]^ The solid 
line is the evolution of the length with time as predicted 
by the numerical solution of equation 
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Figure 9: Time evolution of the average liquid weight on 
each node of a network of crossed fibers making angles 
4.2°, 10.5° and 90°. The sketches illustrate the fiber 
configurations and the initial morphology of the liquid. 
The fiber diameter is 2 a = 0.30 mm, the initial liquid 
(DMTS) volume is 2.7 lJ^^. 


independent. A volume of 2.7 of DMTS is deposited 
on each node. The total amount of liquid placed on the 
network is sufficiently large to measure the mass change 
with time. Thus, the network is placed on a precision 
scale and the weight is recorded over the drying period. 
The experiment is repeated for different angles between 
the fibers by shearing the network. In tuning the angles, 
we modify the liquid morphology at the nodes: column, 
mixed and drop. 

The time evolution of the liquid weight on each node 
is shown in figure The results show that small angles 
enhance significantly the drying rate, which is consistent 
with the results obtained with our model system of two 
crossed fibers since the liquid tends to be in the column 
state for which the evaporation rate is increased. The 
total evaporation time for the mixed conformation and 
the drop are about two and three times the evaporation 
time of the column, respectively. 

Consequently, for applications in which the drying 
time needs to be reduced, we show that a network of 
fibers with small tilt angles performs better. Thus, it 
could be interesting to use fiber networks whose fibers 
have a preferred orientation. Such anisotropy can be, for 
instance, obtained by shearing an anisotropic material 

m- 


4 Conclusion 

In this paper, we present an experimental study on the 
drying of a small liquid volume placed at the node of 
two crossed fibers with different tilt angles. In this con¬ 
figuration, the evaporation of the liquid volume exhibits 
rich dynamics associated with the different morphologies 
of the liquid. Indeed, the drop, column and mixed mor¬ 


phology evaporate at different rates. During the drying 
process, the liquid can explore several of those morpholo¬ 
gies depending upon the initial volume and the angle be¬ 
tween the fibers. Our experiments show that the time of 
drying varies significantly between a column from which 
the liquid evaporates rapidly due to the large surface 
area and the large aspect ratio, and a drop for which the 
small surface area leads to low evaporation rates. To in¬ 
crease the drying rate, it is thus favorable for the liquid 
to be in the column morphology. This morphology can 
be achieved by decreasing the tilt angle between fibers 
as shown here. 

This study provides scaling laws that capture the dry¬ 
ing dynamics of the three morphologies under the as¬ 
sumption of a uniform evaporation flux. Further studies 
would be necessary to investigate more accurately the 
effect of the contact lines on the evaporation flux in this 
particular geometry [3 [3 [To]. We show that the be¬ 
havior documented on two fibers can also apply to more 
complex geometries, such as networks of fibers in which 
nodes are independent. The discussion above focuses on 
drying in a fairly simple geometry, that of two crossed 
fibers. Substrates of higher complexity would render the 
drying dynamics even more complex with possible inter¬ 
actions between neighboring liquid volumes. 
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